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I. INTRODUCTION 

The 60 GHz millimeter wave oxygen absorption line may be utilized to 
observe the temperature of the atmosphere without the problem of water 
vapor attenuation from clouds which would otherwise tend to mask areas of 
the atmosphere from radiometric measurement. To observe wind patterns it 
is desired to measure the temperature of the oxygen with a resolution of 
10 km horizontally as well as vertically. 

A satellite based millimeter wave radiometer can obtain this informa¬ 
tion on more of the atmosphere than a ground based system because of the 
satellite's unobstructed view of the earth. 

A geostationary satellite would provide continuous coverage of a large 
but restricted area of the earth. Angular resolution would be provided by 
the antenna system, and radial resolution would be provided by sweeping the 
frequency of operation across the width of the oxygen absorption line. 

Since the satellite is in a fixed position relative to the earth's 
surface, a continuous communication link with a ground station v?ould permit 
placement of more of the real time control and data processing on the 
earth, thus minimizing satellite payload, providing a potential substantial 
cost savings. 

A polar orbiting satellite would permit tomographic techniques to be 
employed providing coverage of the entire earth, but with only intermittent 
communications with ground stations, and satellite based data processing 
requirements would be substantially greater. In this case single frequency 
operation could be utilized, simplifying hardware requirements since 
tomography would provide the necessary three dimensional resolution. 

Previously, techniques considered to accomplish angular resolution for 
either of the above orbits involved the use a large mechanical structure 


antenna, such as a parabolic dish, utilizing gain to obtain resolution. 



The image of the atmosphere would then be reconstructed utilizing some form 
of beam scanning technique. However, aperture synthesis interferometry, 
long used by radio astronomers, could be utilized instead. This technique 
has the advantage of providing the same resolution of a large dish antenna, 
but instead utilizing an array of much smaller antennas having a maximum 
total spacing equal to the diameter of the dish antenna. Thus a substan¬ 
tial mechanical size and weight advantage may be obtained, resulting in 
less payload cost. 

Aperture synthesis differs from a phased array in that aperture 
synthesis utilizes individual antenna element spacing rather than total 
antenna gain to achieve the same angular resolution. Also aperture synthe¬ 
sis observes the entire field of view at once, whereas a phased array (as 
well as a large parabolic dish antenna) would require some form of beam 
scanning technique. A phased array achieves resolution in the same way as 
a parabolic dish antenna, through the use of high gain, and thus still 
requires the use of a large mechanical structure with a total aperture 
equivalent to that of the dish. 

However, to achieve a substantial payload cost savings over a large 
mechanical structure antenna, an aperture synthesis array will not have the 
same gain because of its much smaller total effective antenna capture area. 
However, signal to receiver noise ratio reduction from reduced gain may be 
corrected, for a slowly varying object to be Imaged radiating thermal noise 
power, by first utilizing as wide a receiver band width as possible, and 
then utilizing integration of the cross-correlated signals. This is the 
technique used by radio astronomers for observation of cosmic thermal 
sources in the "continuum", or broad bandwidth, mode. 

Modern aperture synthesis radiometers require large and complex 
digital processing. However, by placing most of this computer capability 
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on the ground, the payload advantage can be maintained with a communica¬ 
tions link to a ground station. Much of the real time, "synchronous", 
computer control could be combined with the "asynchronous" image mapping 
computer system which would be necessary in any case. Also the error 
compensation flexibility offered by a computerized system may prove to be 
of advantage over an analog beam scanned antenna, particularly in the high 
resolution case. 

Electronic phase stability is also a critical requirement of an 
aperture synthesis instrument. However, techniques can probably be adapted 
from radio astronomy to meet this requirement. Millimeter wave aperture 
synthesis radio telescopes are currently in existence. 

This preliminary report summaries the results of effort performed from 

ylTeS' 

January through January, 1986 toward a study of the application of aperture 

A 

synthesis to this requirement. Chapters II through V describe the 
technical details and requirements of an aperture synthesis interferometer 
in general. Chapter VI uses this theoretical background to calculate the 
approximate antenna requirements for a linear array interferometer in 
geostationary orbit. Chapter VII discusses a possible system concept in 
more detail for a spun linear array aperture synthesis interferometer in 
geostationary orbit. A possible ground-based experiment is also discussed 
which may also prove useful toward an application involving the detection 
of wind shear at airport run ways. Near term future goals of my research 
are also discussed. 

The main emphasis of my research is to proceed toward a system design 
of such a radiometer, and relation of that system design to hardware 
requirements. My previous hardware design experience at the National Radio 
Astronomy Observatory's Very Large Array Radio Telescope Project should 
prove beneficial in this regard. 
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Most of the work performed thus far has been related to gaining an 
understanding of aperture synthesis, ambiguities, and array requirements. 
A preliminary literature search involving more than 60 documents has also 
been conducted. 


A 



II. TWO-ELEMENT INTERFEROMETER WITH COMPLEX CORRELATOR 


A microwave interferometer is an instrument that utilizes the inter¬ 
ference pattern generated by multiplying the signals from two antennas 
together to measure angular direction. This section deals with the geo¬ 
metry of such an instrument, along with the calculation of the output or 
"visibility function" of two classes of interferometer, "simple" and "delay 
tracking". The delay tracking case provides an advantage over the simple 
case in that it minimizes the required field of view in the case of the 
interferometer baseline rotating in the plane of an object to be imaged of 
relatively small total angular dimension. First, the visibility functions 
for the simple and delay tracking cases are calculated in terms of delays 
in Parts A and B. The two visibility functions are then reduced to func¬ 
tions of angle in Part C. 


A. SIMPLE INTERFEROMETER "VISIBILITY" FUNCTION IN TERMS OF DELAYS 

A simple two-element interferometer consisting of two isotropic, 
non-directional, antennas spaced a distance S, measured in meters, is shown 
in Figure II-l. All of the analyses in this paper shall utilize the MKS 
system of units. 

A single received frequency of m in radians/sec shall be assumed for a 
point source at infinity. While the simple point source with a single 
frequency analysis provided here does not provide adequate equations to 
define resolution, bandwidth effects, or image reconstruction, it can be 
used to advantage in understanding field of view and antenna array require¬ 
ments. It is also useful for comparing the field of view performance of a 
simple interferometer to that of a delay tracking interferometer. 

To cover both stationary and rotating situations, the baseline (the 
line connecting the two antennas) is assumed to be rotating relative to the 
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point source so that the angle to the point source from the normal to the 
baseline is given by iJ^Ct) in radians as a function of time. 

The delay in arrival at the second antenna relative to the first 
antenna is given by T(t) in seconds as a function of time. Therefore, if 
the signal present at the first antenna is given by cos mt, where m is the 
received frequency in radians per second, and t is time in seconds, then 
the signal at the second antenna will be given by cos (ji)(t-T(t)) . Note that 
in this simple case, the phase of the received signal at the first antenna 
is assumed to be 0. For the simple two antenna case shown here, this is 
satisfactory, since only relative measurements betv/een the two received 
signals are to be considered. 

Assume a complex correlator is employed as shown in Figure II-l. The 
complex correlator is required to uniquely define the argument of the 
multiplier outputs, since it provides both the sine and cosine of that 
argument, both orthogonal functions. The real and imaginary multiplier 
output calculations shall be handled separately. 

Real Multiplier Output 

Let M , = the output of the real multiplier, 
real 

M , = cos ut cos U)(t-T(t)) 
real 

where t = time in seconds 

uj = received frequency in radians/sec 

T(t) = time delay of arrival of the signal at the second antenna 
relative to the first antenna in seconds. 

From the trigonometric identity, cosa cosB = Jcos(a-S) + icos(a+6), 

M , = Jcos(DT(t) + Jcosa)(2t-T (t)) 
real 

Drop the high frequency, 2u), term with a low pass filter to obtain the 
real part of the complex "visibility" function, 
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= icosa)T(t) . 


Imaginary Multiplier Output 


The Imaginary part of the visibility function may be obtained by 
delaying the phase of one input to the imaginary multiplier by Tr/2 radians, 
Thus, the cos ojt output from the first antenna is delayed in phase by ■ir/2 
radians to provide cos ((i)t-TT/2) . Utilizing the trigonometric identity 
cos(a-T:/2) = sina, 

cos(u)t—ir/2) = sinojt 

Let = the output of the imaginary multiplier. 

M. = cosuCt-T (t)) . sinojt 
imag I' ^uj 

using the trigonometric identity, 


cosasing = Jsin(ct+g) - isin(a-B) 


isina3(2t-T(t)) 
imag __ 


-isina)( 




< ^ V/^ 


ri 


Again, drop the high frequency term with a low pass filter to obtain the 
imaginary part of the complex "visibility" function, 

” -Jsinoj(-T(t)) 

imag 

from the trig identity, 

sin(-ci) = -sina , 

V. = isinu)T(t) 
imag 

Therefore, the complex "visibility" function, expressed in 

terms of delays as: 

V , = V , + jV. 

complex real imag 

Complex = + jisin,,x(t) 

If the — delay were Introduced into the second antenna's signal path 
instead, it can be shown that the conjugate of the above function will be 


obtained. 
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B. DELAY TRACKING INTERFEROMETER "VISIBILITY" FUNCTION IN TERMS OF DELAYS 
A delay tracking two element interferometer consisting of two iso¬ 
tropic (nondirectional) antennas spaced a distance s measured in meters is 
shown in figure II-2. Note the addition of the delay tracking device to 
the simple interferometer described previously. The significance of this 
device will be explained in section III. ^ ^ 

A single received frequency of 'jt radlans/sec shall be assumed for a 
point source at infinity. 

To cover both stationary and rotating situations, the baseline (the 

line connecting the two antennas) is assumed to be rotating relative to the 

point source so that the angle to the point source from the normal to the 

baseline is given by ijj(t) in radians as a function of time. 

The delay at the second antenna relative to the first antenna is given 

by xCt) in seconds as a function of time. Therefore if the signal present 

at the first antenna is given by cosmt, where o) is the received frequency 

in radians per second, and t is time in seconds, then the signal at the 

second antenna will be given by cosmCt-x(t)). 

The first antenna's output, cosmt, is delayed in time by in 

seconds by the addition of the delay tracking device. Therefore the signal 

to the left side of the first multiplier is given by coscoCt-x^^^Ct)). Note 

that X called the reference delay, is also a function of time 

ret 

because of the relative motion of the baseline. Again assume a complex 
correlator is employed as shown in figure II-2. The calculation of the 
real and imaginary multiplier outputs shall again be handled separately. 


Real Multiplier Output 


Let the output of the real multiplier 

M = cosaj(t-x ^(t)) • cosm(t-x(t)) 

ir63iX ir0X 
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where t = time in seconds 

u = received frequency in radians/seconds 

T(t) = time delay of arrival at the second antenna relative to the 

first antenna in seconds 

T^^^(t) = reference time delay in seconds 

From the trigonometric identity, 

cosacosB = icos(a-6) + icos(a+0), 

M T = Jcosa)(T(t) - T ,(t)) 
real ref 

+ icosa)(2t - T ^(t) - T(t)) 
ret 

Drop the high frequency, 2to, term with a low pass filter to obtain the 
real part of the complex "visibility" function, 

Imaginary Multiplier Output 

The imaginary part of the visibility function may be obtained by 
delaying the phase of the input to the imaginary multiplier by ^ radians. 
Thus the cosa)(t - output from the delay tracking network is 

delayed in phase by radians to provide cos [tD(t-T^^^(t)) - . Utilizing 

the trigonometric identity, 
cos(a- = sina, 

cos[a)(t - T^g^(t)) - |] = sina)(t - 

Let ~ output of the imaginary multiplier 

M. = cosa)(t-T(t)) • sina)(t-T c(b)) 

imag ref 

Using the trigonometric identity, 

cosasin3 = isin(a+8) - Jsin(a-3) 

M = Jsin(jo(2t-T(t)-T (t))-isina)(T (t)-T(t)) 

imag ref ret 

Again, drop the high frequency term with a low pass filter to obtain the 
imaginary part of the "visibility" function, . 
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''l.ag ■ 

From the trigonometric identity, 

sin(-a) = -sina, 

V. = Jsina)(T(t) - t ,.(t)) 
imag ref 

Therefore the complex "visibility" function, expressed in terms 

of delays as; 

V , = V , + jv. 

complex real imag 

V ^ = Jcosu)(T(t) - T _(t)) 

complex ref 

+ jisinai(T(t) - ■^^gf(t)) 

Again if the ^ delay were introduced into the second antenna's signal path 
instead, it can be shown that the conjugate of the above function will be 
obtained. 


C. VISIBILITY FUNCTIONS IN TERMS OF ANGLE AND TIME 

Referring to figure II-3 for angle relationships to time delays, the 
visibility function can be expressed in terms of time and angle instead of 
delays. 

Let il;(t) in radians be the angle measured from the normal to the 

baseline to any point source as a function of time. 

Let Ip ^(t) in radians be the angle from the normal to the baseline to 
ref 

a reference point in the center of the object to be mapped (imaged), again 
as a function of time. It is not necessary that the reference point be 
exactly in the center of an object to be mapped, but if it is, a minimum 
relative field of view will result. "Field of view" is defined as the 
maximum angular dimension of the source, and is further discussed in 
chapter III. 

Let Alp in radians be the offset angle from the reference angle 

ip ..(t). Note that it is not a function of time, 
ref 
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= 'j'j.gf(t) + Ai|) 


where tjj(t) = angle in radians to any point source as a function of time 
= reference angle in radians as a function of time 
All) = offset angle in radians of point source relative to refer¬ 
ence direction 

Let a)^g£ in radians per second be the angular rate of rotation of the 
interferometer baseline relative to the reference direction. Therefore, 

ill ir(t) “(1) 

^ref ref 

where reference angle in radians as a function of time 

u) ^ = angular rate of rotation in radians/second, 

ref 

t = time in seconds 
combining the two equations, 
l|/(t) = t+Aii) 

where i|;(t) = angle in radians to any point source 

u ^ = angular rate of rotation in radians/second 
ref 

t = time in seconds 


Ai|) = offset angle in radians 

Let s in meters be the spacing between the two antennas. Let c be the 
speed of light in meters/sec. Therefore the time delay along the baseline 
is given by s/c in seconds. 

Let in seconds be the relative time of arrival delay of the 

signal at the second antenna relative to the first antenna from an imagin¬ 


ary point source at the reference angle From trigonometry, then 

“ 7 sin(m ,*t) 
ref c ref 

where T^g£(t) = relative delay in seconds from an imaginary reference 


point source 


s = spacing in meters 


c = speed of light in meters/sec 
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^ = rate of rotation of baseline relative to reference 

ref 

direction in radians/sec. 
t = time in seconds 

Let T(t) in seconds be the time of arrival of the signal at the second 

antenna relative to the first antenna from any point source. 

Similarly then, 

T(t) = - sln(c» - • t + AiJ;) 
c ref 

where T(t) « relative delay in seconds from any point source 
s = spacing in meters 
c = speed of light in meters/sec 

w , = rate of rotation of baseline relative to reference direction 
ref 

in radians/sec 
t = time in seconds 

Aij' = offset angle in radians to any point source 
The above expressions for delay may now be substituted into the 
previous visibility equations expressed as functions of delay to obtain the 
visibility functions in terms of rate of rotation, ^rid offset angle, 

A<|>. 


Simple Interferometer 

V , = icos(*)T(t) + iJsin<«)T(t) 

complex 

substituting for t(t), 

qOJ 

V T = Jcos[— sin(<i> ^*t+A>^)] 

complex c ref 

e(l) 

+ jj|sin[— sin(t^ ^•t+A'l')] 
c ref 

Let f in Hertz be the received signal frequency corresponding to ui in 
radians/second. 
a)=2Ttf 

Substituting for 
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V . = icos[^^^ sin(a) ,*t+Ai|;)] 

complex c ref 

For convenience, express the visibility function in terms of the baseline 

spacing in wavelengths, s^, which is dimensionless. Let X in meters be the 

wavelength of the received frequency f, in Hertz, 
s 

" X 

where = baseline spacing in wavelengths, dimensionless 
s = baseline spacing in meters 

X = wavelength of the received frequency, in meters. 

But ^ ~ ^ 
sf 

s = — 

X c 

where s, = baseline spacing in wavelengths, dimensionless 
A 

s = baseline spacing in meters 
f = frequency of source in Hertz 
c = speed of light in meters/sec. 

Substituting back into the previous visibility equation yields: 

V , = icos[2iTs, sln(aj ^‘t+Ail))] 

complex X ref 


+ jJsin[2iTs^ sin(•t+Ai);) ] 


where V . = the complex visibility function for the simple interfero- 

complex 

meter 

s = baseline spacing in wavelengths 
A 

03 = rate of rotation in radians/sec 

ref 

t = time in seconds 
Aijj = offset angle in radians 


Delay Tracking Interferometer 

V = icososCrCt) - T -(t)) 

complex ret 

+ j isin o)( T(t) - 
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substituting for T(t) = ^ sin( ’t+Aij;) 

and T *(t) = — sin(u ^'t), 
ref c ref 

V , = Jcos(— [sin((u £ •t+A'I') 

complex c ref 

+ jjsin(— [sin(t»)^^^’t+Ail;) 
Utilizing the trigonometric identity 


sin(a)ref ’t) ]) 
sin(u)^^^ -t) ]) 


sinQ - sin3 = 2cosi (“+3)sini » 

V = Jcos[^^ cos((»i ^'t+J Aip)sin(i AiJ/) ] 

complex c ref 

+ jjsin['^^ cos (*t+i Aij^) sin (J AiJj) ] 

Again substituting u = Zirf, 

V = ^cos cos((i) ^’t+i Ai|;) sin(i A'i') ] 

complex c ref 

+ jjsin[^—^ cos((i)^^^’t+J Aii))sln( J Aij)) ] 

Again substituting 

V T = icos[Ars, cos(<») ^'t+J A4 j) sin( J AiJ;) ] 

complex A ref 

+ jJsin[4Trs, cos (to ’t+i Aij;) sin(i Ai|;) ] 

” rer 

where ~ complex visibility function for the delay tracking 

interferometer 


Sj^ = baseline spacing in wavelengths 

to - = rate of rotation in radians/sec 
ref 

t = time in seconds 
At() = offset angle in radians 

It can be shown that the delay tracking visibility function reduces to the 

simple interferometer visibility function for either t=0, the reference 

direction coincident with the normal to the baseline, or to ^=0, the 

ref 

stationary case. 

Let to ‘t-O, then 
ref 

V , = Jcos [Alts-.cos(i Aij;)sin(i AiJj) ] 

complex ^ 

+ j |sin t Alts i^cos ( i Aii»)sin(i AiJ)) ] 


lA 



Using the trigonometric identity cosasina = j!sin2a, 

V , = Jcos [2Trs .sinAV] + jj sin [2 its sin AH'] 

complex X A 
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III. SPATIAL APERTURE SYNTHESIS 


A. SPATIAL FREQUENCY 

The synthesis of a large aperture antenna to achieve high angular 
resolution imaging in the plane of the baseline of the interferometer can 
be accomplished by varying the spacing of a two element interferometer and 
measuring the visibility function as a function of spacing. Spatial 
frequency is then used with the Nyquist sampling criteria to determine the 
number of antennas in a "spatial" aperture synthesis interferometer 
required to meet a given field of view and resolution requirement. The 
visibility functions for the simple and delay tracking interferometers are 
given by: 



A:j; = offset angle from reference angle in radians 

If the spacing in number of wavelengths, s^, is varied while holding 

all other variables constant, the real and imaginary components of the 

visibility functions will vary sinusoidally with s^ with a radian spatial 

frequency of m .. In other words, the visibility functions of both cases 
S A 

of interferometer may be expressed as, 

V = Jcos((Dg^ • Sj^) + jjsin(u)gj^ ■ s^) 
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Therefore the radian spatial frequency, w for both cases of interfero- 

S A 

meter may be expressed as: 

Simple Interferometer 

u) , = 2Trsin(a) ^'t + Ai|)) 
s A ref 

Delay Tracking Interferometer 

*^sA ~ ^■fcos(io^^^‘t + i Ai|)) sin(i Alp) 

The spatial frequency, f^^ may be defined as 
0 ) 

, s A 


sA 2-n 

where f , is the spatial frequency, dimensionless 
s A 

<i) is the radian spatial frequency in radians 

S A 

2IT has units of radians. 

Hence, the visibility functions of both cases of interferometer may be 
expressed as, 

V = Jcos(2irf ,*s,) + jisin(2TTf ,*Si) 

S A A S A A 

where f . = spatial frequency, dimensionless 
S A 

Sj^ = number of wavelengths spacing, dimensionless 

2fT = constant with dimension of radians 

Therefore the spatial frequency for both cases of Interferometer may be 

expressed as: 

Simple Interferometer 

f , = sln( 0 ) ,'t + Alp) 

s A ref 

Delay Tracking Interferometer 

f 1 =2cos(w , *t + J Ai|)) sin (i Alp) 
s A ret 

The spatia]. frequency, fg is both a function of time and offset 

angle Aip, assuming the rotation rate is a known constant. Since time, 

t, and rotation rate known quantities, the spatial frequency can 

be measured to determine the position of a source, the offset angle, Aip. 
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This can be accomplished by either moving one antenna relative to the 
other to vary s^^, or by properly spacing many elemental antennas such that 
the sinusoidal visibility function is sampled adequately. In the latter 
case, the Nyquist sampling rate which determines the minimum number of 
antennas required is dependent on the maximum spatial frequency encoun¬ 
tered. A comparison of maximum spatial frequencies for both the simple and 
delay tracking interferometer cases is presented in figure III-l. in 

radians is defined as the total maximum angular extent of the object. The 
peak angular field of view, in radians is defined as one half the 

maximum angular extent of the thermal source to be imaged. Thus it repre¬ 
sents the full range of Ai|; values of sources that the interferometer may 
observe. For small angular extent sources in a cold sky background the 
source itself may limit the maximum angular extent encountered. Let in 
radians be defined as the total maximum angular extent of the object to be 
imaged. In general, should be limited by individual elemental antenna 

gain, such that other sources (e.g. interference signals from other satel¬ 
lites, cosmic radio sources, the sun, etc.) are prevented from reaching the 
correlator. 

The relationship of to the angular geometry is shown in figure 

III-2. Both the "restricted" case, for observation of sources of "re¬ 
stricted" field of view and the "general" case, for sources at all possible 
angles are pictured. Note that the general case is valid for all angles no 
matter what possible ambiguities may exist. An example of a practical 
system which exists that validates this approach to the angular extent of 
the "general" case is that of a narrow aperture direction finder where two 
interferometers are placed perpendicular to each other to provide cosine 
and sine components of direction. This type of system is unambiguous, but 
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MAXIMUM SPATIAL FREQUENCIES 


SIMPLE f 

s\ 

= sin(u) -•t+Aib) 
ret 



CASE 

U) iT* t 

ref 

Aijj 

^max 

^sX max 

General 

Rotating 

Any 

Any 

1 

Restricted 

Rotating 

Any 

^^fv 

1 

General 

Stationary 

0 

Any 

1 

Restricted 

Stationary 

0 

^’*'fv 

sin(Ai|»f^) 

DELAY TRACKING f 

sX 

= 2cos(aj^^^. t+iAt(i)sin(iAip) 


CASE 

“ref- ^ 

A)l/ 

^max 

f 

sXmax 

General 

Rotating 

Any 

Any 

2 

Restricted 

Rotating 

Any 

^’<'fv 

2sin(iAij# 

General 

Stationary 

0 

Any 

1 

Restricted 

Stationary 

0 


sin(Atp 


FIGURE III-l 
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is capable of observing signals at all angles, even though by synnnetry each 
interferometer is ambiguous. Ambiguities are discussed in more detail in 
Chapter IV. 

The importance of delay Cracking is indicated in figure III-l by a 

comparison of f^^ max between the simple and delay tracking interferometers 

for the restricted rotating case. For sources of small angular extent in a 

rotating system, delay tracking is seen to provide an enormous advantage 

over the simple type of interferometer, reducing the spatial frequency by a 

factor of 2 sin(iAi|( ). 

r V 

B. RESOLUTION AND ARRAY LENGTH 

From the resolution required for a given application, the maximum 
linear dimension, or aperture, can be calculated for an aperture synthesis 
interferometer. 

Resolution is defined as the ability to make distinguishable the 

H 

individual parts of an object. Christiansen and foghorn state that the 

Rayleigh convention for the resolution of an optical instrument is not 

1 

suitable for a radio telescopef Johnson and Jasik state that two sources 

will be resolved if the half power beamwidth of the antenna system is less 

than one half the beamwidth between the first nulls.^ This is stated to 

usually be the case. Therefore the half power beamwidth may be utilized as 

1 

the limiting resolution case of a radio telescopef Verschuur and Kellerman 
give the minimum half power beamwidth of a minimum redundancy linear array 
(as proposed in section IIIC to minimize hardware requirements) as the 
wavelength divided by the maximum linear dimension of the array.^ 

Therefore, 

^ . = - 

resmin s 

max 
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where \jj = minimum resolution angle in radians 

resmin ° 

X = wavelength in meters 

s = maximum linear dimension of the array in meters 
max ■' 

In terms of the maximum array spacing in number of wavelengths, 
s 

_ max 
^ ^ax X 

where = maximum linear dimension of array in numbers of wavelengths, 

dimensionless 

s = maximum linear dimension of array in meters 
max •' 

X = wavelength in meters 

Therefore, the requirement for the maximum dimension of the array is given 

ty. 

^ I 
Xmax i(i 

resmin 

where “ maximum linear dimension of array in numbers of wavelengths, 

dimensionless 


. = minimum resolution angle in radians 

resmin 

and the maximum array length or aperture can be calculated using. 


s 

max 


X 

’j' 

resmin 


where s maximum array length in meters 

max ® 

X = wavelength in meters 

. = minimum resolution angle in meters, 

resmin 

Rogers states that angular resolution varies with direction by the 

component of the maximum baseline normal to the source direction.^ Figure 

III-3 shows that the component to the baseline is given by 

s, = s, cos(m ,-*t) 

Amaxc Xmax ref 
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where s= component of maximum antenna spacing normal to the direction 
of the source in numbers of wavelengths, dimensionless 
s^ajj “ maximum antenna separation in numbers of wavelengths, 
dimensionless. 

= angular rate of rotation of array in radians/sec. 
t = time in seconds. 

, =—L_ 

res s, 

Amaxc 


where 4/ = resolution angle in radians in any reference direction. 

^^maxc ~ component of maximum antenna spacing normal to the direction 
of the source in numbers of wavelengths, dimensionless. 

Note that the resolution angle has been selected to be in the reference 
direction usually to the center of the object to be imaged. For a 
satellite mapping the Earth's atmosphere this is sometimes called the 
"well" resolution angle. Therefore, 


res 


S. ‘COSCU x: *t) 
Amax ref 


where “ resolution angle in radians in any reference direction. 

^Amax ~ maximum antenna spacing in numbers of wavelengths, 
dimensionless. 

(i) ^ = rate of rotation in radians/second of antenna array in plane 

ref 

of source. 

t = time in seconds. 

Note that the minimum resolution angle occurs at ui ^'t = 0, therefore for 

° ref 

the calculation of array requirements it shall be assumed that either the 
normal to the antenna baseline points toward the source or that at some 
time this situation still occurs, in order to minimize microwave as well as 
digital processing hardware requirements. 
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C. LINEAR ARRAY ANTENNA REQUIREMENTS 

Utilizing information on the field of view and resolution required 
along with the previous information on the maximum spatial frequency to be 
encountered one can predict the number of antennas required for a linear 
array aperture synthesis interferometer. 

If a finite number of elemental antennas are used in a linear array, 
the spatial domain will be sampled as shown in figure III-A which shows a 
zero redundancy configuration and spacings for A antennas. It is assumed 
that the output from each antenna is multiplied by the outputs from all of 
the other antennas to obtain the spacings required. A plot of the real 
part of the sinusoidal visibility function is superimposed on the diagram 
to illustrate the Nyquist sampling requirement in relation to the spatial 
frequency. The Nyquist sampling theorem requires that each visibility 
function spatial period must be sampled at least two times to insure no 
loss of information. 

Note that in practice, the required sampling rate will be somewhat 

higher. Adapting Stremler's time signal analysis to the spatial signal 

analysis here, two limitations raise the required sampling rate.^ 

First the spatial signal that is being sampled is spatially bounded, 

from s, . to s. , not from to +", and a spatially limited signal can 
Amin Amax 

never be spatially frequency limited. Second a perfect elemental antenna 

cannot be built that would be capable of totally rejecting spatial 

frequencies above f , 

s Amax 

Both of these limitations will thus require a somewhat higher sampling 

rate to achieve given desired worst case aliasing errors. Therefore the 

following analysis is in reality a best case analysis for the minimum 

number of required samples, as well as for A , , the minimum number of 

min 

required elemental antennas. 
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The Nyquist sampling theorem can be expressed albegraically as: 


2(Jmax < p 

N . s^min 

min 


where 


’Amax 


N 


rain 


= spatial width of each sample, dimensionless 


s, = maximum antenna spacing, aperture, in numbers of 
Amax 

wavelengths, dimensionless 

N . = minimum number of non-redundant antenna spacings, 

min 

dimensionless 

^s^min ~ "'inimum spatial period of visibility function, dimension¬ 
less 


2s, 


Therefore N . > 


'max 


min P 


sAmin 


since f 


1 


sAmax P 


sAmin 


N > 2 s f 
min Amax sAmax 

where N . = minimum number of non-redundant antenna spacings, dimension- 

min 

less 

Sa = maximum antenna spacing, aperture, in numbers of wavelengths 
''max 

dimensionless 

f = maximum spatial frequency of visibility function, dimension 


less 


The maximum number of "non-redundant" spacings or minimum number of 
samples to satisfy the sampling theorem for 1 to 4 antennas is given by 
A. T. Moffett® as 

N - MAzii 
min 2 
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where N ^ = maximum number of non-redundant spacings available from a 

min 

given number of antennas or minimum number of samples required 
to satisfy the Nyquist sampling theorem, dimensionless. 

A = total number of elemental antennas. 

Unfortunately, 1 to 4 antennas are the only linear arrays having no 
redundancy^. Above four antennas, the spacings obtained will not be evenly 
spaced, no matter what algorithm is chosen to space them. However "minimum 
redundancy" arrays can be chosen. An example of a minimum redundancy array 
of 5 antennas is shown in Figure III-5. In this case one of the baseline 
spacings is redundant thus only providing 9 nonredundant spacings, whereas 
10 total possible spacings are available from A(A-l)/2. 

Therefore, a redundancy number, R, is defined using the equation, 

N = 
min 2R 

where = minimum number of non-redundant spacings required to satisfy 

the Nyquist sampling criteria. 

A = number antennas 
R = redundancy 

In the above example of five antennas, R is calculated to be 10/9 or 1.11. 

Figure III-6 tabulates values of R for minimum redundancy.® For 

larger linear arrays of antennas than 11, the redundancy has been proven to 

be no less than 1.217 and no greater than 1.332®. 

Note that the equality for will also be valid for non-integer 

spacings if one takes N . to be the ratio of the greatest spacing, the 

min 

aperture, to the greatest sized incremental spacing. The greatest sized 
incremental spacing is the worst case spacing to meet the given Nyquist 
sampling requirement for the maximum visibility function spatial frequency 
requirement based on the angular field of view of the source to be imaged. 
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Utilizing some algebra and the quadratic formula, an equation can be 
found for the number of antennas, A, to produce a required minimum number 
of non-redundant spacings, to satisfy the sampling theorem require¬ 

ment . 

N = 
min 2R 

2RN . = - A 

min 

A^ - A - 2RN , =0 

min 

Using the quadratic formula to solve the above equation, 

a(A)^ + b(A) + c = 0 , A = ^ ^ 

2a 

a = 1, b = -1, c = -2RN ^ 

min 


1 ± / 1 - 4(1){-2RN . ) 



1 ± 'i 1 + 8RN ^ 

. min 

A = -2- 

Since only positive and integer values of A are of interest, the minus sign 
can be dropped, and the integer function can be used to determine A . , the 
minimum integer number of antennas required to produce at least a required 
minimum number of spacings. The integer function will round down, 
therefore 1 must be added to it to insure an adequate number of antennas. 
Therefore, 

1 + / 1 + 8RN . 

A . = 1+INT[-r-2^] 

min 2 

where A . = minimum integer number of antennas (dimensionless) required 

min 

to produce at least a required minimum number of non- 

redundant spacings, N . . 

min 
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INT = integer function 
R = redundancy, dimensionless 

N^in = required minimum number of spacings to satisfy the Nyquist 
sampling theorem requirement. 

For a large number of antennas, there will almost always be a greater 

number of non-redundant spacings available than are required because of the 

sum of 1 with the integer function as well as the square root. Previously, 

it was found that N . > 2s ■, f . . Therefore, it is valid to drop the 

min Amax sAmax ' 

greater than sign for all practical purposes for large arrays. Therefore 
let, 

N . = 2s, f , 

min Amax s Amax 

for substitution into the previous antenna formula. Substituting for N . 

min 

then yields. 


A . = 1+INT [ 

min 


1 + / 1 + 16Rs 


Amax s Amax 


] 


where “ minimum integer number of antennas, dimensionless, required 

to adequately sample a given spatial frequency 
INT = integer function 

R = redundancy, dimensionless, see Figure III-6 for values 
^Amax ^ dimension of array in number of wavelengths, 

dimensionless 

^sAmax ° maximum spatial frequency, dimensionless 

However, the resolution angle requirement determines the 

maximum linear dimension or aperture of the array. This relation was 

previously derived as 

^ I 
^ Amax 'I' 

resmin 
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Substituting into the minimum antenna equation. 


X ^^x 

A , = 1+INT [ -5— 

min 


f is dependent on the field of view requirements as stated in Figure 

s Amsx 

TII-l. Thus the dependence of A . on Aiii, may be derived using whichever 

min fv 

one of the eight cases listed in the table of figure III-l is appropriate. 
Because of the expense of microwave and digital processing hardware and 
ambiguity considerations presented in section IV, only three of the eight 
cases would be considered practical. The first two of these are included 
in the restricted field of view stationary case, which results in the same 
f for either the simple or the delay tracking interferometers where, 

S Ain^x 

f , = sinAi()- 

s Xmax f V 

The third practical case is the restricted field of view rotating case of 

the delay tracking interferometer, which would be the minimum spatial 

frequency case for a rotating interferometer. In this case, 

f = 2siniAiJ;- 

s Xmax fv 

Substitution results in the final equations relating the number of antennas 
required for given worst case peak field of view and resolution angle 
requirements. Note that these should be considered best case design 
equations for the reasons stated earlier. 


Restricted Stationary 

Simple or Delay Tracking Interferometer 


A . = 1+INT [ 

min 



16Rsin A')'£^ 


resmin 


1 
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Delay Tracking Interferometer 



where A . = minimum integer number of antennas, dimensionless, for 

min 

given field of view and resolution requirements 
INT = integer function 

R = redundancy of linear array, dimensionless, see figure III-6 
for values 


= worst case peak field of view in radians 

il; . = minimum resolution angle in radians, 

resmin 

These two equations predicting the number of antennas required for a linear 
array, along with the equation relating the maximum linear dimension or 
aperture of the array. 


s 

max 


X 

'I' 

resmin 


where s = maximum spacing or aperture of array in meters 


X = wavelength in meters 


tL . = minimum resolution angle in radians 

resmin 

./III I ■ (yt'l) 

define the maximum complexity and size of the array. Possible practical 
situations are presented in Chapter VI. 


D. ELEMENTAL ANTENNA REQUIREMENTS 

The beamwidth of each elemental antenna will generally be approximate¬ 
ly equal to the field of view, since this will result in the minimum 
spatial frequency visibility function to be sampled. Elemental antenna 
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gain will affect the signal to noise ratio received, but in general, this 
is not significant, if the time required for integration can be tolerated, 
since integration of the cross-correlated signals may be employed to 
generally obtain any desired signal to noise ratio for imaging of thermal 
noise sources. 

Thus an aperture synthesis interferometer can achieve a substantial 
mechanical structure advantage over a parabolic or phased array antenna for 
the same resolution requirement. 

In a practical situation severe constraints may be placed on the 
elemental antennas in terms of how rapidly they reject interfering signals 
outside of the field of view. These signals may arise from the sun, other 
satellites, cosmic radio sources, etc. Since they are outside the field of 
view they will not be sampled adequately in the spatial frequency- domain 
and may be aliased into the image. The elemental antenna requirements thus 
have to be determined by an error analysis done on the actual application 
and type of interferometer employed. 

E. FOURIER TRANSFORM PROCESSING 

While the previous analyses covered only a single frequency point 
source, in a practical situation infinitely many such sources are present 
simultaneously. A discrete Fourier transform is then utilized to sort out 
the different resolution elements from which an image of the entire object 
being observed is constructed. This provides the basis for the "spatial" 
approach to aperture synthesis. Rotation of a two element interferometer 
with time may also be utilized but this has not been examined as yet. 
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IV. AMBIGUITIES 


Ambiguities in the angular position of a source measured with respect 
to the simple two element interferometer may result from several indepen¬ 
dent causes. Ambiguities due to geometry, the type of correlator utilized 
with a given baseline spacing, undersampling for a given field of view, and 
Fourier transform "windowing" are discussed in this section. 


A. GEOMETRY 

Ambiguities in the delay function T(t) for the total angle cj^^^-t+Atp 
exist due to the symmetry of the geometry of the simple two element inter¬ 
ferometer with isotropic (non-directional) antennas. Figure IV-I illus¬ 
trates the ambiguities due to geometry for both the three dimensional case 
as well as the more restricted two dimensional case. The analyses for both 
cases assume that the time of arrival difference between the two antennas 
can be measured unambiguously by the rest of the interferometer system. 

This assumption isolates the ambiguities due to geometry from ambiguities 
due to other causes. Figure IV-1 is drawn consistent with standard 
engineering notation for a right hand spherical coordinate system, r is 
defined as the radius to the source in meters, 0 is defined as the angle 
from the positive z-axis to the source, and (fi is defined as the angle from 
the X axis to the projection of the source in the x-y plane. 

In this analysis, both the desired source and any undeslred ambiguous 
sources are assumed to be at infinity. The coordinate system used in the 
previous planar analyses is maintained with the use of defined as the 
angle measured from the perpendicular to the antenna baseline to the source 
in the plane defined by the two antennas and the source. The delay function 
T is defined as the delay in the time of arrival of the signal from the 
source at antenna B relative to antenna A. 
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1. Three Dimensional Ambiguities 

From Figure IV-1, the desired source is located at coordinates (r=“, 

&=6 , <J)=(|) ). The delay function t is given by: 
s s 

T = — cos 6 
c s 

where t = delay of time of arrival of signal at antenna B relative to 
antenna A 

s =■ spacing between the two antennas in meters, 
c = speed of light in meters/sec. 

0^= angle to source measured from the positive z axis. 

Note that t is independent of (j). Therefore, any signal at coordinates 

(r=% 0=0 , ) is ambiguous with the desired source at (r=», 0=0 ,4)=(j) ) 

s s s s 

since exactly the same time of arrival delay r is achieved. Directions of 

sources resulting in Identical delays in three dimensions are indicated by 

the conical surface containing antenna B and the circle shown in Figure 

IV-1. Note that these directions are the only ambiguities resulting from 

geometry. Sources at r=“ and contained within only the one conical surface 

given above will result in Identical delays. Any source at (r=“, 0 =tt- 0 ,any 

s 

(Ji) will result in a negative delay and so will not be ambiguous. 

The ambiguous sources described by coordinates (r=“>, 0= 0^, (()^(|)^) result 
from the cyclindrical symmetry of the simple two-element array. Thus the 
addition of more elements to the linear array does not aid this problem 
since array symmetry is maintained. Also directional antennas cannot 
eliminate three dimensional ambiguities since can be arbitrarily close 

to (b=(b , the direction of the desired source. 

Two methods, however, may be implemented to limit three dimensional 
ambiguities to the two dimensional case. 

In one method, a planar array of antennas may be utilized instead of a 
simple linear array. A second method would be to rotate the linear array 
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in a plane not including the source and preferably perpendicular to the 
source. In both cases cylindrical symmetry would be eliminated and only 
the two dimensional ambiguities described in the following section would 
exist. 

It could be said that the National Radio Astronomy Observatory "Very 
Large Array" aperture synthesis radio telescope utilized both techniques. 
Its planar "Y" configuration was effectively rotated with respect to 
astronomical sources by the Earth's motion due to the latitude of the VLA 
site in New Mexico. If the VLA site had been located at the north or south 
poles, only a simple linear array would have been required. However, at 
the lower latitude in New Mexico, only limited rotation was achieved for 
non-polar sources and thus three linear arrays in a "Y" shaped configura¬ 
tion provided near complete baseline coverage of radio sources in a range 
of directions. 


For satellite based observations of the earth from a geostationary 
satellite, the rotational method would probably be the most economical, 
since a single linear array could be used to synthesize a planar array, and 


generate all baselines required to meet resolution requirement^^-In—fact, 
a-near irrfl-nite number-of baselines eeuld be ac hi e v ed*at no additional cost 


V /V 

in microv7ave hardware and=w©«^d also thus minimiz^ satellite payload. 


2. Two Dimensional Ambiguities 

In this case, ambiguities are restricted to the plane containing the 

two antennas and the source. Thus is a constant defining that plane. 

Therefore, the coordinates of the desired source are given by 

(r=", 6=9 ,(|)=(t) ), where <t> is considered a constant. The delay function t 
s s s 

is again given by: 
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g 

T = — COS 6 

c s 

where t = delay of time of arrival of signal at antenna B relative to 
antenna A. 

s = spacing between the two antennas in meters, 
c = speed of light in meters/sec 

angle to source measured from the positive z axis. 

From the trigonometric identity 
cos(oi) = cos (-a) 

g 

then, T = — cos(-0 ) 
c s 

Therefore (r=", 6=-6 ) is ambiguous with (r=“, ft=0 ) since the same 

s s s s 

delay results from the two directions. 

Note that the three dimensional case Includes the two dimensional 

case. In three dimensions, 

(r=", 0=0 , ) is ambiguous with (r=«“, &=0 , ) 

s s s s 

therefore 

(r=», 0= 0^, 4)=(j)^+iT) is ambiguous with (r=“, 6=0^, (|>=()>g) 

but 

(r=“, 0=0 ,4>=4>+Tf) is equivalent to (r”®®, 0=-0 , )• 

s s s s 

Therefore (r=“, 0=-0 , 4>=<|> ) is ambiguous with (r=", 0=0 ,ij)=<|) ), the 

s s s s 

previously derived two dimensional case. 

g s 

Also, note that the delay t = —cos0 is equivalent to t = —slni|j pre¬ 
sented in chapter II. 

s a 
T = —COS 6 
c 

S / ^ I \ 

T = -COS(^lJ<) 

c z 

Using the trigonometric identity 
cos(^a) = sinci, 
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g 

therefore, t = —simj;. The ambiguities derived here in terms of 6, can 
similarly be converted to terms of if desired. 

B. CORRELATOR AMBIGUITIES 

In this section the complex correlator presented in Chapter II is 
shown to not measure time delay uniquely for a continuous source at any 
angle of arrival. However the maximum antenna spacing in wavelengths 
can be calculated for unambiguous operation of two element simple or delay 
tracking interferometers. The problem of ambiguous operation for large 
spacing Interferometers is then shown to be utilized to advantage in 
unambiguous operation of aperture synthesis interferometers. 

The complex correlator visibility functions for a two element inter¬ 
ferometer in terms of time delays were previously derived in Chapter II as 

Simple Two-Element Interferometer 

V - = J cos(OT(t) + jJsinu)T(t) 

complex 

where V , = complex visibility function 

complex 

(ti = received frequency in radians/seconds 
T(t) = time delay of arrival at the second antenna relative to 
the first antenna in seconds 
Delay Tracking Two-Element Interferometer 

V = JcostD(T ^(t)-T(t)) + j |sinu)( T ,(t)-T(t)) 

complex ref ref 

where V , = complex visibility function 

complex 

0 ) = received frequency in radians/seconds 
T(t) = time delay of arrival at the second antenna relative to 
the first antenna in seconds 

Note that in both cases, the cosine and sine functions are orthogonal 
functions which are capable of uniquely defining the arguments only over a 
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TT radian range individually or a Zir range when used together. Thus if 
|a)T(t)|<Tr then the complex correlator will provide an unambiguous visibil¬ 
ity function for the simple interferometer. 

Delay tracking can be seen to provide an advantage in this limitation 
if the peak field of view is small, since x ^(t)WT(t). 

If lw(T^^^-r Ct))| <it, 

then the complex correlator will provide an unambiguous visibility function 
for the delay tracking interferometer. If the peak field of view of the 
object is small then T^^j(t)~T(t) and a much larger x (t) can be accommo¬ 
dated before an ambiguity will result. 

The effect of antenna to antenna spacing can most easily be examined 
in general by writing the complex visibility function as a function of 
spacing and spatial frequency. In this case. 


V = |cos(2Trf ,s ) + jisln(2Trf s ) 

complex sX X sX X 

where V , = complex visibility function 

complex 

f^^ = spatial frequency, dimensionless 
s = spacing between the two antennas in number of wavelengths. 


dimensionless 


2ir = constant with dimensions of radians 
The cosine and sine functions are orthogonal functions which are 
capable of uniquely defining the argument, 2Tif^^s^ , only over a 2t^ radian 
range when used together. Therefore | 2Tr f ^^s^ | <-ir will result in unambiguous 
operation. Therefore, the range of antenna spaclngs allowed for an unam¬ 
biguous visibility function is given by 


1 



By utilizing the table of maximum spatial frequencies derived in 
Figure IIl-l, along with the formula for unambiguous spacings, a table of 
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maximum spacings for each case of two-element interferometer and possible 
field of view can be derived. This table of unambiguous spacings is 
presented in Figure IV-2. 

The delay tracking interferometer is seen to have an advantage over 
the simple Interferometer in the restricted rotating case for objects with 


a small field of view, The maximum spacing in this case may be 

increased by 2| sin(lAi|; —Tf from that of the simple interferometer and yet 
still maintain unambiguous operation. 


However, the maximum spacing of an aperture synthesis interferometer 


must be set by the desired resolution angle, ijj . , and not the much 

resmin 

larger peak field of view angle, As described in section III-B, the 

maximum spacing required for an aperture synthesis interferometer is: 



but |i) , for a high resolution instrument. Therefore, 

^res ^fv 

s, 

Xmax 

for any of the cases shown in Figure IV-2. This will certainly be ambigu¬ 
ous in a two element fixed spacing interferometer, and will result in many 
cycles of the sinusoidal visibility function as shown in Figure IV-3. 

While the correlator will produce an ambiguous visibility function 
output for a two-element interferometer, this problem is utilized to 
advantage in an aperture synthesis interferometer. By varying the spacing, 
through either physically moving one antenna relative to the other, rotat¬ 
ing the two antennas in the plane of the source to change the projection of 
the baseline perpendicular to the directory of the source, utilizing 
discrete elemental antennas at strategic spacings to evenly sample the 
visibility function, or a combination of the above, the visibility function 
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UNAMBIGUOUS ANTENNA SPACING 
FOR 2-ELEMENT INTERFEROMETER 


SIMPLE 

f , = sin(u) t-hAi|)) 

sX ref 



CASE 

w t 

ref 

A\|i 

max 


General 

Rotating 

Any 

Any 

< 1 

Restricted 

Rotating 

Any 

^^fv 

< h 

General 

Stationary 

0 

Any 

< J 

Restricted 

Stationary 

0 

"^fv 

2] sin(A\!^^^) 1 


DELAY TRACKING f , = 2cos(w t+iAif;) sin (iAiJ;) 
- sX ref 


CASE 

ID t 

ref 

Aij; 

max 


General 

Rotating 
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Any 

< i 

Restricted 

Rotating 

Any 

^-^fv 
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A] sin(jAi|)^^) 1 


General 0 Any < \ 

Stationary 


Restricted 

Stationary 


0 



2| sin(Aii»^^) I 


FIGURE IV-2 


















will appear as a periodic sinusoidal function versus spacing. By utilizing 
a Fourier processor complex objects may be mapped without correlator 
ambiguities in an aperture synthesis interferometer. 


C. UNDERSAMPLING FOR A GIVEN FIELD OF VIEW 

As discussed in section Ill-C, the Nyquist sampling theorem requires a 
band limited periodic spatial signal unbounded in the spatial domain. Both 
of these conditions do not hold for the reasons previously presented. Thus 
an error analysis is required for a given system to determine the actual 
number of samples and elemental antenna pattern required to limit the worst 
case error due to sidelobes resulting from aliasing. 


D. FOURIER TRANSFORM V7IND0WING AND MAPPING TECHNIQUES 

The inverse Fourier transform required for processing the spatial 

domain is limited by the spaclngs to Techniques are available 

called "windowing" which minimize sidelobe effects due to this limitation 

on the inverse Fourier transform and must be investigated, again with a 

suitable error analysis. It may be possible to adapt mapping techniques 

.C 

such as "plean" from radio astronomy to aid in this process. 
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